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Multi- wavelength light curves of bright gamma-ray blazars (e.g., 3C 454.3) reveal strong correla- 
tions across wavebands, yet striking dissimilarities in the details. This conundrum can be explained 
if the variable flux and polarization result from both (1) modulation in the magnetic field and rela- 
tivistic electron content imparted at the jet input and (2) turbulence in the flow. In the Turbulent 
Extreme Multi-Zone (TEMZ) model being developed by the author, much of the optical and high- 
energy radiation in a blazar is emitted near the 43 GHz core of the jet as seen in VLBA images, 
parsecs from the central engine, as indicated by observations of a number of blazars. The model 
creates simulated light curves through numerical calculations that approximate the behavior of tur- 
bulent plasma — modulated by random fluctuations of the jet flow — crossing a cone-shaped standing 
shock system that compresses the plasma and accelerates electrons to highly relativistic energies. A 
standing shock oriented transverse to the jet axis (Mach disk) at the vertex of the conical shock can 
create a variable nonthermal seed photon field that is highly blueshifted in the frame of the faster 
jet plasma, leading to highly luminous, rapidly variable 7-ray emission. 



I. INTRODUCTION 

The variability of blazars from radio to 7-ray fre- 
quencies has long been a challenging problem that 
theoretical models have yet to solve. The causes of 
this futility include complexity in the data, which in 
the past have been undersampled in both time and fre- 
quency. The routine all-sky survey mode of the Fermi 
Large Area Telescope (LAT) , plus the concerted effort 
by many research groups to provide commensurate 
data at other wavebands, are combining to provide 
new data sets that are sufficiently rich to both test 
existing scenarios for blazar emission and stimulate 
the design of new models. The hope is that finding 
the correct explanation for the variations will allow us 
to use the data to infer the structure of the jet as close 
to the black hole as possible. 

The author is leading a collaboration dedicated to 
amassing extensive multi-waveband fluxes, polariza- 
tion, and parsec-scale images over time. This has led 
to the recognition that the current array of models 
with, at most, a few emission zones are unable to re- 
produce some prominent aspects of the observational 
data. This work summarizes the relevant character- 
istics of the data and introduces a new numerical, 
many-zone scheme that can potentially explain many 
aspects of the observed behavior of blazars. 



II. LOCATION OF GAMMA-RAY FLARES 

Successful modeling of blazar emission requires that 
we identify the location(s) where the flares occur. In 
order to do this, the author and his collorators use 
their extensive multi-waveband data to determine the 
timing between when a high-energy flare occurs and 
when a superluminal knot passes through the "core" 
on millimeter-wavelength Very Long Baseline Array 



(VLBA) images. (The core is an approximately sta- 
tionary, bright, compact feature at one end of the 
jet on sub-milliarcsecond (parsec) scales. At mil- 
limeter wavelengths, it has observed properties that 
are similar to those expected from a standing conical 
shock. 0,Q) They employed this method in an analysis 
of BL Lac in 2005 to determine that there are multiple 
sites of X-ray and optical flares in the jet: upstream of 
the core and within the core.[lO( It appears that flares 
occur as enhanced plasma (higher magnetic field and 
particle density) in the jet flow passes through these 
regions. Evidence for this is provided by the quasar 
PKS 1510—089, which exhibited a series of 7-ray and 
optical flares in early 2009 as the optical polarization 
position angle x rotated systematically by ~ 720°. As 
this rotation ended, a bright superluminal knot passed 
through the core at (to within the uncertainties) the 
date of an extremely bright optical/7-ray flare. The 
last flare therefore took place either in or very close 
to the core. 

Further evidence that the majority of flares at 7-ray 
(and usually other) frequencies originate in or near 
the millimeter-wave core comes from the coincidence 
in timing of millimeter- wave and 7-ray flares. Q How- 
ever, the 7-ray peak usually occurs during the rise of 
the mm-wave flux, leaving the possibility that the 7- 
ray event takes place somewhat upstream of the core. 
The author's group monitors the structure of the jets 
of 35 Fermi-detected blazars with the VLBA at 7 mm. 
After three years of this program, preliminary analysis 
finds that 43 7-ray flares were simultaneous within er- 
rors with the appearance of a new superluminal knot 
or a major outburst in the core at 7 mm. This com- 
pares with only 13 cases in which either a 7-ray flare 
occurred without a corresponding mm-wave event, or 
a mm-wave flare or superluminal knot ejection had no 
7-ray counterpart. (Four of these blazars had other 7- 
ray flares that did correspond to a mm-wave event.) 
In addition, five blazars that were quiescent at 7-ray 
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energies were also quiescent at 7 mm. 

The evidence is therefore clear that most (but not 
all) of 7-ray flares occur near or in the mm-wave core. 
The challenge is then to create a model within this 
context that is capable of reproducing the observed 
characteristics of the multi-waveband variability of 
blazars. 



III. OBSERVED PROPERTIES OF BLAZAR 
VARIABILITY 

Over the past decade, a number of detailed studies 
of long-term, well-sampled light curves of blazars at 
various regions of the electromagnetic spectrum have 
uncovered some patterns that should be incorporated 
in theoretical models. One is the power spectral den- 
sity (PSD) of the flux variations, which in blazars is 
generally has a red noise spectrum, meaning a power 
law with steep slope. [l|, U [l6[ This suggests that a 
random process modulates the emission. For syn- 
chrotron and inverse Compton radiation, the random 
process must govern fluctuations in the magnetic field 
and density of relativistic electrons. Turbulence is a 
likely cause of such fluctuations. This hypothesis is 
supported by the degree of linear polarization, which 
indicates the level of disorder of the magnetic field 
vectors across the emission region. If we approxi- 
mate the turbulence in terms of N discrete cells, each 
with its own magnetic field orientation, the mean de- 
gree of polarization should be p ~ p m ax-^ - > where 
Pmax ~ 0.75 is the fractional polarization in the pres- 
ence of a uniform magnetic field. [2] A typical optical 
polarization is of order 10%, hence we can surmise 
that there are ~ 60 turbulent cells present at optical 
wavelengths. The polarization tends to be lower at 
millimeter wavelengths. Q This can be explained if the 
emission region is larger than at optical wavelengths, 
encompassing more turbulent cells. Finally, the time- 
scale of variability is generally shorter at higher fre- 
quencies between the radio and optical bands. [15] This 
can result from averaging of the fluctuations over more 
cells at longer wavelengths. 

Close examination of multi-waveband outbursts 
provides further evidence for chaotic behavior super- 
posed on more systematic trends. The quasar 3C 279 
is a prime example. Figure Q] presents 2.4 years of light 
curves at 7-ray, X-ray, and optical wavelengths, and 
in the 7 mm core (as seen in VLB A images). Times 
when a new superluminal knot passes the 7 mm core 
coincide with various events: (1) the beginning of the 
largest 7-ray and optical outburst, (2) an "orphan" X- 
ray flare, and (3) a flare that occurred simultaneously 
at 7-ray, optical, and X-ray energies. In addition, a 
major outburst in the core started before the 7-ray 
flare at RJD 5460-5480 (where RJD = JD-2450000). 
Furthermore, a period of quiescence occurred contem- 
poraneously at 7-ray, X-ray, and optical bands be- 
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FIG. 1: Multiwavelength light curves of the quasar 3C 279 
from August 2008, when regular 7-ray monitoring by 
Fermi began, to the end of 2010. Note the overall sim- 
ilarity of the 7-ray and optical light curves but also the 
differences in the detailed variations. The vertical arrows 
in the bottom panel indicate times when new superlumi- 
nal radio knots passed the core on 7 mm VLBA images. 
The flux of the core plotted includes the blended flux of 
the new knot until it separates from the core by more than 
0.1 milliarcsec. Data are from Marscher et al. (in prepa- 
ration). 



tween RJD 5300 and 5440. It is clear that the rela- 
tionships among the wavebands are complex. Other 
blazars with multiple outbursts — e.g.,PKS 1510—089 

in 2009 IHG3 and 3C 454 - 3 in late 2010 (Wehrle et 
al., in preparation) exhibit similar complexities along- 
side obvious (and statistically significant) correlations. 

The difference between 7-ray and optical fluctua- 
tions is particularly noteworthy. For magnetic field 
strengths between ~ 0.2 and 1 G, as is usually inferred 
in the optical synchrotron-emitting regions, electrons 
of the same energies (thousands times the rest-mass 
energy) radiate at both optical wavelengths and 0.1- 
10 GeV energies. The implication is that something 
besides the number of such electrons is changing to 
produce the variations in flux. The synchrotron flux 
depends on the magnetic field strength B as well as 
its direction, number of radiating electrons V rc , and 
Doppler factor 5, while the inverse Compton flux de- 
pends on only the latter two. If the seed photons 
that are scattered arise from outside the jet, they are 
blueshifted in the plasma frame, in which case the 7- 
ray emission is more sensitive to S. Since the 7-ray 
flux usually (but not always) fluctuates with higher 
amplitude than does the optical flux, either B is anti- 
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FIG. 2: Multi-waveband flux and degree P and electric- 
vector position angle \ of optical polarization vs. time 
for BL Lac in 2011. Note rapid variations in P and x> 
indicative of a turbulent magnetic field. Data are from 
Marscher et al. (in preparation). 




FIG. 3: Sketch of geometry adopted in the TEMZ model. 
Top: View down the jet, whose cross-section is divided 
into many cylindrical cells. Actual calculations displayed 
in Figs. 4-6 use 271 cells across the jet. Bottom: Side-on 
view. Cylindrical cells appear as rectangles. The conical 
shock compresses the flow & accelerates electrons. The 
Mach disk, if present, is at the axis, oriented transverse to 
the flow. 



correlated with N (which seems unlikely on physi- 
cal grounds), the highest-energy electrons are created 
mainly in cells with magnetic fields pointing almost 
toward the observer (after correction for relativistic 
aberration), or the Doppler factor varies considerably 
across the emission region. The last of these possibil- 
ities is a feature of "jet-in-jet" models involving mag- 
netic reconnection events [6j or relativistically turbu- 
lent motions. [13 These have been proposed to explain 
extremely short time-scales of variability of the TeV 
7-ray emission of some blazars. 

BL Lac (Fig. [2]) provides another example of an 
object whose multi-waveband light curves, polariza- 
tion, and changes in sub-milliarcsecond-scale struc- 
ture support the scenario in which turbulence plays a 
major role in the nonthermal emission. A number of 
7-ray flares are seen, but their relationship to optical 
and X-ray variations is complex. In 2011, the blazar 
became paticularly active at all wavebands starting 
about RJD = 5700, as a new superluminal knot was 
passing through the core at 43 GHz (as occurred dur- 
ing all of the main 7-ray flares seen in Fig. [5]). The 
optical polarization became highly chaotic until just 
after a rapid TeV flare was observed by VERITAS. 18] 
Another 7-ray flare near RJD = 5865 occurred during 
a major X-ray and radio outburst after another su- 
perluminal knot appeared, while the optical flux was 
fluctuating rapidly and x(°pt) was again chaotic. Tur- 
bulence is a natural way for the jet to produce these 
types of variations. [14 



IV. TURBULENT EXTREME MULTI-ZONE 
(TEMZ) MODEL 

The author is developing a model in which much 
of the optical and high-energy radiation in a blazar 
is emitted parsecs from the central engine. The main 
physical features are a turbulent ambient jet plasma 
that passes through a conical standing shock wave 
(recollimation shock) in the jet. The model can gen- 
erate short time-scales of optical and gamma-ray vari- 
ability by restricting the highest-energy electrons ra- 
diating at these frequencies to a small fraction of the 
turbulent cells, perhaps those with a particular ori- 
entation of the magnetic field relative to the shock 
front (Summerlin & Baring 2012). Because of this, 
as well as radiative energy losses as the plasma ad- 
verts beyond the shock front, the volume filling factor 
at high frequencies is relatively low. Such a model is 
consistent with the (1) red- noise power spectra of flux 
variations in blazars, (2) shorter time-scales of vari- 
ability of flux and polarization at higher frequencies, 
(3) mean polarization levels as well as fractional devia- 
tions from the mean that are higher at optical than at 
lower frequencies, (4) apparent rotations in polariza- 
tion position angle that are really just random walks 
of the projected magnetic field direction, (5) breaks in 
the synchrotron spectrum by more than the radiative 
loss value of 0.5, and (6) flares that are often sharply 
peaked or contain multiple maxima, both of which 
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FIG. 4: Two simulated spectral energy distributions 
(SEDs) from the TEMZ code. In both cases, the left side 
of the SED is from synchrotron radiation, while the right 
component is from inverse Compton scattering. For the 
latter, the seed photons are from IR from dust (dashed 
curve) and synchrotron plus SSC emission from the Mach 
disk (solid curve). 



are not reproduced by single- or few-zoned models. 
The dependence of items 2-4 on frequency is directly 
related to the change in spectral index beyond the 
break, according to the model. 

The model includes synchrotron radiation, inverse 
Compton scattering (IC) of seed photons from hot 
dust [observed in 4C21.35Q, and IC of synchrotron 
plus synchrotron sclf-Compton (SSC) radiation from 
relatively slowly moving plasma in a Mach disk (shock 
oriented perpendicular to the axis in cylindrically 
symmetric jets) at the conical shock's vertex. If 
the Mach disk is present, this is the dominant — and 
variable — source of seed photons, since the emission 
is Doppler boosted in the frame of the turbulent cells. 
The combined effects of non-uniform electron energy 
distribution, different magnetic field orientations for 
different turbulent cells, and light-travel delays often 
cause time lags and/or single waveband ("orphan") 
flares (see Figs. 4-6 for sample output from prelim- 
inary computations). Thus, the model can qualita- 
tively reproduce some of the puzzling features of ob- 
served multi-waveband light curves of blazars. In ad- 
dition, the discrete 7-ray/optical correlation function 
is similar to that observed in blazars, with zero mean 
lag but considerable spread about this value. 

In the future, the author plans to develop the 
TEMZ model by adding the following capabilities: 




Time steps 



FIG. 5: Sample TEMZ light curves: optical synchrotron 
radiation (black, dashed curve) and inverse Compton 7- 
ray flux. Seed photons for the latter are synchrotron and 
first-order SSC emission from the Mach disk. Red arrow 
points to an orphan 7-ray flare. 



1. "Peculiar" velocities of the cells relative to the 
mean flow, as expected for turbulence or magnetic 
reconnections. This is probably necessary to repro- 
duce intra-day 7-ray variability with high apparent 
luminosities, [a, [l7| 

2. Full SSC calculation, with seed photon field for 
each cell including retarded-time emission from all 
other cells. This is relevant mainly to BL Lac objects 
with 7-ray to synchrotron integrated flux ratios < 1. 
It will require translation of the code from Fortran to 
C++ and parallelization to run on a supercomputer. 

3. Refine the calculation of pair-production absorp- 
tion at high 7-ray energies and synchrotron self- 
absorption at low frequencies. The code currently 
does this crudely in order to save computation time. 
These sources of opacity affect the spectral energy dis- 
tribution (SED) at millimeter wavelengths and at 7- 
ray energies above a few GeV. 

4. Split each cell into ~ 10 sub-cells for a more refined 
calculation, to allow intra-day variability. The current 
scheme assigns uniform physical properties across each 
cell. Because of the acute angle of the conical shock, 
this requires long, thin cells (see Fig. 3). 

5. Calculate the polarization of X-rays from inverse 
Compton scattering, in the anticipation that this can 
be measured by GEMS (launch scheduled for 2014) 
for some blazars during 7-ray/X-ray flares to add ad- 
ditional constraints on models. 
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FIG. 6: Sample TEMZ light curves: optical synchrotron 
radiation (black) and inverse Compton 7-ray emission 
(red). Seed photons for the latter are from thermal in- 
frared emission by a dusty torus. 

V. CONCLUSIONS 

While analyzing light curves as a series of distinct 
flares can be a convenient way to study the physical 
conditions that exist at different flux states of a blazar, 
the statistics of light curves indicates that noise pro- 
cesses dominate. The degree and variability of linear 
polarization also lead directly to a picture in which the 
magnetic field has a major chaotic component, which 
can result from turbulence. The higher polarization 
and stronger variability at shorter wavelengths further 
suggests that the volume filling factor of the emission 
decreases with frequency. While this can arise from 
radiative losses after electrons are energized at a shock 
front [HI, [lij , the break in the synchrotron spectrum 
is usually greater than the change from spectral in- 
dex a to (a + 0.5) expected from this process. A 
sharper break can occur if electrons can be acceler- 
ated to higher energies only over a smaller fraction of 



the total volume. 

The author is incorporating these requirements into 
a qualitatively new paradigm for both the multi- 
waveband variability of blazars and the appearance 
of superluminal radio knots. Stochastic variations in 
the jet flow — changes with time of the magnetic field 
and electron density at the jet input, perhaps driven 
by fluctuations in the accretion flow or magnetic field 
configuration in the central engine — combined with 
turbulence within the flow cause the variability in 
emission as the plasma crosses standing shocks on par- 
sec scales in the jet. Seed photons for inverse Comp- 
ton scattering can come from thermal infrared emis- 
sion by a hot dust torus, synchrotron radiation from 
the turbulent plasma, or synchrotron plus first-order 
SSC radiation emitted beyond a Mach disk, which is a 
small, strong shock oriented perpendicular to the jet 
axis that slows the flow down to subsonic speeds. The 
emission from the Mach disk is highly blueshifted in 
the frame of the plasma that crosses the conical shock 
that extends across most of the jet cross-section. 

The rapid variability in the TEMZ model occurs be- 
cause of the small volume filling factor of the highest- 
frequency emission, as well as the fluctuations in mag- 
netic field, electron density, and maximum electron 
energy from one turbulent cell to another. The addi- 
tion of velocity structure in the turbulence will accen- 
tuate the variability on the shortest time-scales. After 
developing the numerical TEMZ code further, the au- 
thor will run it over a wide range of parameter space 
in order to compare the simulated light curves and dy- 
namic SEDs with those that are observed in blazars. 



Acknowledgments 

This research is supported in part by NASA 
through Fermi grants NNX08AV65G, NNX10AO59G, 
and NNX11AQ03G, and by NSF grant AST-0907893. 
The VLBA is an instrument of the National Radio 
Astronomy Observatory, a facility of the National Sci- 
ence Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



[1] A. A. Abdo et al, Astrophys. J. 722, 520 (2010). 
[2] B. J. Burn, Mon. Not. Royal Astron. Soc. 133, 67 
(1966). 

[3] T. V. Cawthorne, Mon. Not. Royal Astron. Soc. 367, 
851 (2006). 

[4] R. Chatterjee et al, Astrophys. J. 689, 79 (2008). 
[5] F. D. D'Arcangelo et al, Astrophys. J. Let. 659, L107 
(2007). 



[6] D. Giannios, D. A. Uzdensky, and M. C. Begelman, 

Mon. Not. Royal Astron. Soc. 395, L29 (2009). 
[7] S. G. Jorstad et al, Astronom. J. 134, 799 (2007). 
[8] J. Leon-Tavares, E. Valtaoja, M. Tornikoski, A. Lh- 

teenmki, A., and E. Nieppola. Astron. Astrophys. 

532, A146 (2011). 
[9] M. P. Malmrose, A. P. Marscher, S. G. Jorstad, 

R. Nikutta and M. Elitzur, Astrophys. J. 732, 116 



eConf C110509 



G 



2011 Fermi Symposium, Roma, May 9-12 



(2011). 

[10] A. P. Marscher et al, Nature 452, 966 (2008). 
[11] A. P. Marscher et at, Astrophys. J. Let. 710, L126 
(2010). 

[12] A. P. Marscher et at, in 2009 Fermi Symposium 
(eConference C091122 ), ed. W . N. Johnson and D. 
J. Thompson. larXiv:1002.0806V l (2010). 

[13] A. P. Marscher and W. K. Gear, Astrophys. J. 298, 
114 (1985). 

[14] A. P. Marscher, W. G. Gear, and J. P. Travis, in Vari- 
ability of Blazars ed. E. Valtaoja and M. Valtonen 
(Cambridge U. Press), 85 (1992). 

[15] A. P. Marscher and S. G. Jorstad, in Fermi Meets 



Jansky: AGN in Radio and Gamma Rays, ed. T. 
Savolainen, E. Ros, R. W. Porcas, and J. A. Zen- 
sus (Max-Planck-Institut fur Radioastronomie, Bonn, 
2010), p. 171. 

[16] I. M. McHardy, X-Ray variability of AGN and rela- 
tionship to Galactic black hole binary systems, in The 
Jet Paradigm, ed. T. Belloni, Lecture Notes in Physics 
794, 203 (2010). 

[17] R. Narayan and T. Piran, Mon. Not. Royal Astron. 
Soc. 420, 604 (2012). 

[18] R. A. Ong, ATel #3459 (2011). 



eConf C110509 



